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2. Stability of the Analog Low-Cost Large Space Telescope Due to 
Quantization 

2-1 . Introduction 

The objective of this chapter is to conduct an investigation on 
the pointing stability of the low-cost Large Space Telescope (LST) 

System under the Influence of quantization at various locations of 
the system. Only analog model of the LST system is considered 
in this chapter. ie stability of the digital LST system with quan- 
tization is considere In Chapters 5 and 6. 

The block diagram of the simplified analog LST system with 
quantizers is shown in Fig. 2-1. 

The describing function method is used to determine the condition 
of self-sustained oscillation In the LST system due to the effect of 
each of the three quantizers acting alone. The present analysis 
considers only one quantizer at a time. 

Let the analog describing functions of the quantizers Qp, Q^, and 
Q R be represented by Np, N^, and N R , respectively. In general, the 
describing function of a quantizer nonlinearity is a function of the input 
amplitude, E; quantization level, h; and the number of quantization 
levels, n, which depends on E. 

For the LST system shown in Fig. 2-1, the "characteristic equations" 
of the system when each one of the quantizers is acting, are given as 


follows: 


Displacement Quantizer Q q : 


V 3 ♦ N D (K p s ♦ Kj) + K r s 2 - 0 (2-1) 

Torque Quantizer Q^.: 


J y s 3 ♦ N t (K r s 2 ♦ K p s ♦ Kj) ■ 0 (2-2) 

Rate Quantizer Q r : 


J y s 3 + K p s + Kj ♦ N^s 2 • 0 (2-3) 


These equations can be conditioned by dividing both sides of the 
equations by the terms that do not contain the describing function, so 
that th-i stability r uation is expressed in the form of 

1 + NG(s) - 0 (2-4) 

where N is the describing function, and G(s) is a linear transfer 
function. The condition of self-sustained oscillations is found by 
investigating the oossible intersections between the trajectories of 
-1/N and G(s) in the complex p'ane. 

For the equations given in Eqs. (2-1), (2-2) and (2-3), the equiva- 
lent transfer functions are 


ty‘> 


V + K i 

S 2 (J„S ♦ K r ) 


(2-5) 


t 








Consider that a quantizer has the Input-output relation as shown 
In Fig. 2-2. Let the Input to the quantizer be a sine wave, 

e ( t ) ■ E s 1 n wt 

The quantization level is h, and let the magnitude of E be such that 

(2n - Dh . , E < .Ua^Jin 

where n is a positive integer. A typical output of the quantizer is 
shown in Fig. 2-3. 

The fundamental component of the Fourier series expansion of y ( t ) 
is 

y 1 (t) = Y 1 sin wt 


where 


Y, -i 

I 71 


tt/2 


y ( t ) s i n wt dwt 


Evaluating the last integral, we have 


Y i = 2- 

1 77 


^n 4 r* 72 

y(t)sinwtdwt + - 1 nhsinwtdwt 


a-. 


71 J 


a. 


n-1 


■I .1 


fO 


i = l 


i+1 


ih sin wt dwt + 


If 


tt/2 


nh sin wt dwt 


Ol: 




quantizer corresponding to a sinusoidal input. 


h ■ - T '["‘“HI • co$ “t3 * — cos “n 

Since 



Expanding the right-hand side of Eq. (2-13), we can show that 


(2-18) 


(2-19) 


( 2 - 20 ) 





Figure 2-4 shows the plot of N(E/h) as a function of E/h. We 
shall show In the following that 


Lan N( E/h ) - 1 

f- 


Letting x ■ E/h, Eq. (2-20) Is written 

, 2 ) 1/2 


>“ ■ k j, l' • M } 


Expanding the quantity Inside the summation sign In the last 
equation, we have 


Or 


N(x) 


-4- ? 1 

"X fa [ 

, in 

n 2 ifi n 4 

i 

roj 
x | 

• 

2xj 8(x ’ 2xJ 


3 (1 

1 l 6 15 fl l 


‘ 48 [ X * 

• 384[x ‘ 2 




N(x) 'nh 


l 

2 ? 


1 - 


ll 


1 


• • 4^' ‘ 'I 


15 

384 X * 




( 2 - 21 ) 


( 2 - 22 ) 


(2-23) 


(2-24) 


Taking limit as x approaches infinity and n approaches infinity, in 
the last equation, we get 
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Um N(x) • Um £- [ 

y-*oo X"** 9 * 1*1 

n-*" n-« 



— ? I 

2x Z 8 X 4 
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n 



1 

7 


s1n'\ 


1 

o 


- 1 (2-25) 

Since the describing function of the quantizer nonlinearity is 
always a real number, the function -l/N(E/h) in the magnitude (db) 
versus phase coordinates will lie on the -180-degree axis for all 
values of E/h. The plot in Fig. 2-4 shows that the magnitude of 
-l/N(E/h) is infinite for 0 < E/h < 0.5. For 0.5 < E/h < 0.707, the 
plot of -l/N(E/h) extends from infinity to -2.09 db along the -180-degree 
axis. Over the range of 0.707 < E/h <_ 1.5, -l/N(E/h) extends from 
-2.09 db to 1.9 db along the -180-degree axis, etc. As E/h approaches 
infinity, the plot of -l/N(F/h) is reduced to the zero-db point. 

Figure 2-5 shows the plot of N(E/h) in magnitude versus E/h and shows 
the multivalued property of the function. 





the quantizer nonli 














Figure 2-6 shows that since the plot of Gp ( s ) does not intersect 
the -l/N(E/h) trajectory in the finite domain, the quantizer 0 Q will 
not cause any self-sustained oscillations in the LST system. 

Both the curves for G T (s) and G R (s) intersect the -180-degree axis 
and the -l/N(E/h) trajectory at 25 db. 

Thus, 


20 ^g 10 |l/N(E/h)| = 25 db (2-26) 

which gives 

N( E/h ) * 0.05012 (2-27) 

Figure 2-4 shows that for this value of N(E/h), 

£ * 0.5 (2-28) 

or 


E s 0.5 h 


(2-29) 
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This means that the quantisers Qy and Q R each Independently may 
cause self-sustained oscillations to occur In the 1ST system. Further- 
more, the self-sustained oscillations will have an amplitude equal to 
one-half the level of quantization h. Another point of Interest Is 
that the Intersect between Gy(s), G R (s), and -l/N(E/h) corresponds to 
n ■ 1 , so that tne quantizers Qy and Q R are essentially actinq as a 
simple relay with dead zone. 

Figure 2-6 shows that the frequency of the oscillations caused by 
Qy Is 1.5 rad/sec, and that caused by Q R Is 6.3 rad/sec. 

It is Interesting to point out that the analysis in Chapter 1 shows 
that the displacement quantizer Q D produces far qreater quantization 
error than the torque and rate quantizers. However, the describing 
function annysis shows that Q D does not cause self-sustained oscillations, 
whereas Qy and Q R may excite oscillations with amplitudes equal to one- 
half the level of quantization, h/2. All these factors must be taken 
into consideration when selecting the quantization level. 
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3. Stability of the Analog Low-Cost Large Space Telescope with 
Reaction Wheel Friction Nonlinearity 

3-1 . Introduction 

The objective of this chapter is to conduct an investigation on 
the pointing stability of the analog model of the low-cost Large Space 
Telescope (LST) System with the reaction wheel frictional nonlinearity. 

Although the LST system is digital, the stability study on the 
analog model will establish a limiting case when the sampling period 
becomes very small, thus providing a check on the results of the digital 
system. 

The block diagram of the digital low-cost LST system, including the 
reaction wheel dynamics, is shown in Fig. 1-1. It has been shown in 
Section 1-1 that the dynamics of the LST system can be simplified. 

Figure 3-1 shows the block diagram of the simplified low-cost LST 
system with the reaction wneel nonlinearity. 

The A of the system shown in Fig. 3-1 is 

A = 1 + K r G 4 + + NG 6 G ? (3-1) 


where N denotes the continuous-data describing function of the reaction 
wheel frictional nonlinearity. 

Substitution of the expressions of G^ , G^, G^, G^, and G^ into 
Eq. (3-1 ) , we have 


K R K I 1 

A- 1 (Kn+T-)- 1 ? 

V p s J s Z 


+ N 


J RW S 


(3-2) 








Setting A to zero, and rearranging Eq. (3-2), we qet 

VrH * 3 * J R«V 2 ♦ * K I» * N V ' 0 

which Is the "characteristic equation" of the system. 
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3-2. Condition of Self-Sustained Oscillations in the Analog LST System 
with Reaction Wheel 

It has been shown [1] that the nonlinear frictional characteristics 
of a reaction wheel can be described by the Dahl model. Therefore, the 
analog describing function of the CMG frictional lonlinearity derived 
in [ 2 ] cai e directly utilized. 

The equivalent transfer function that N sees Is determined from 
Eq. (3-3) by dividing both sides of the equation by the terms that do 
not contain N. We have 

NJ s 

1 + 5 ** 0 (3-4) 

W + K R 5 * V * V 


Thus , 



(s) 


J RW<V' 


V 

T 


+ k r s 


+ V ♦ Kj) 


(3-5) 


The condition of self-sustained oscillation Is investigated by 
plotting G^ju)) and -1/N in the magnitude (db) versus phase coordinates. 
Figure 3-2 shows the plot of G eq (j(i)) of the LST system with 

K r * 3.71 x 10 5 
Kp = 1.65 x 10 6 
Kj » 7.33 x 10 5 
J v * 41822 


0.2 


60 



Figure 3-2. Ggg(juj) and -1/N plots of analog LST with reaction wheel and 
frictional nonlinearity. 


and the plots of -1/N with 


(1) Y “ 8470 (11) Y » 847C0 (111) Y ■ 8470 

T 6F o - 0.0424 T 6F0 - 0.0424 Tq F0 * 0.424 

Since the G eq (jaj) plot never enters the region bounded by -180° 
and -270° In which the plot of -1/N lies, the analog 1ST system will 
not have self-sustained oscillations, due to the reaction wheel 
frictional noni ir^arlty. 
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4. Stability of the Digital Low-Cost Larqe Space Telescope W<*h 
Reaction Wheel Friction Nonlinearity 


4- 1 . Introduction 

In this chapter we shall Investigate the condition of self- 
sustained oscillations of the digital low-cost LST system with reaction 
wheel friction nonl Inearlty. 

The block diagram of the simplified LST system is shown in Fig. 

4-1. The equations written for the outputs of the samplers are: 

♦ # U) - -G A (z)4> e (z) + N(z)G b (z)O rw (z) (4-1) 

0 RW (z) ■ G c (z)<t e (z) - N(z)G d (z)O rw (z) (4-2) 

where N(z) denotes the discrete describing function of the reaction 
wheel friction nonlinearity. Since the friction nonlinearity of the 
reaction wheel can be represented by the Dahl model, N(z) Is identic. 1 
to the discrete describing funct on of the CMG nonlinearity [2], 

The transfer functions In Eqs. (4-1) and (4-2) are 


G A< Z > ** 


G ho G l G 2 G 3 1 




G B (z) = & 


G ho G 2 G 3 
1 + k r g 2 


(4-3) 


(4-4) 











?3 


G 0 (Z) 


.(W 

ir ^2 


The signal flow graph representing Eqs. (4-1) and (4-2) Is shown 
In Fig. 4-2. The A of the system Is obtained from Fig. 4-2, 

A - 1 ♦ G a (z) ♦ N(z)[G a (z)G 0 (z) - G 0 (z) G c (z) ♦ G d (z)] 


Setting A to zero, the last equation can be written In the form 


1 ♦ N(z)G eq (z) - 0 


<W'> 


G a (z)G 0 U) - g b (z)g c (z) * g d (z) 


The Individual transfer functions are evaluated as follows: 


A V s‘(s ♦ a) s 3 (s ♦ a), 


■ W*> * r l 1 ^-^ 

p B J v 2a (z - 1)‘ 


T . 1 z - 1 

a 2 (z - 1) a 3 a 3 (z - e' aT ) 


a “ K R /J v 




1 - e 




4-2. Self-Custalned Oscillations In the Digital 1ST System With 
Reaction Wheel Nonlinearity 


Figure 4-3 shows the plots of G (z) In Eq. (4-9) for various 

eq 

values of n with T as a parameter. The Integer n and sampling period T 
are related to the frequency of oscillation by the following equation: 



(4- 


The following system parameters are used: 

K * 1.65 x 10 6 K 0 - 3.71 x 10 5 

P R 

Kj = 7.33 x 10 5 J y = 41822 

J RW " °’ 2 


Ir Fig. 4-3 the curve for n * 2 extends up to approximately 30 db 
at T = 0.7 sec. Also, as n is increased the curves for G gq (z) approach 
the curve for G eq (s) presented in Fig. 3-2. The plots of G eq ( z ) 
together with the plots of -1/N(z) allow the study of self-sustained 
oscillations in the digital LST system. 

Figures 4-4 through 4-6 show the plots for -1/N(z) for various n. 
In these plots y = 84700 and T^ » 0.424. In all cases, the magnitude 
of the lowest point of the -1/N(z) curve as E approaches 0 is given by 


Lim 

E-*0 


TOT 



( 4 - 


With the given parameters this point is at -83.66 db. The curves 
for n = 2 are shown in F.j. 4-4 and consist of two straight lines at 
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-180° (for 0 <_ ♦ < n/2) and -360° (for */2 < <}> £ n). The plots for 
n ■ 3 and n ■ 4 arc shown In Figs. 4-5 and 4-6, respectively. Several 
values of <J> are plotted In each case to Illustrate the effect of the 
phase of the Input signal. It should be noted that for odd n the curves 

repeat every 180/n degrees starting from <J> * 0, and for even n the curves 

repeat ever: 360/n degrees starting from &, ■ 0. As the Input amplitude E 

goes to Infinity, the curves also go to Infinity and span a region 180/n 
degrees or 360/n degrees wide for odd or even n, respectively. This 
region is centered about the -270° line. Thus, as n goes to Infinity 
these curves approach the -1/N curve of the continuous system shown in 
Fig. 3-2. 

For stability analysis it Is sufficient *o consider only the bounds 
of the -1/N(z) plot for a fixed n. Self-sustained oscillations can 
occur If the G(z) curve corresponding to the same n intersects with 
the F(z) plot. Figures 4-3 through 4-6 show that self-sustained oscillations 
can readily occur in this system with the choice of T determining the 
possible n values which can exist. As the critical regions for higher 
values than n = 4 can be easily visualized. It is apparent that this 
system will have self-sustained oscillations for Integral valcas beyond 
n = 4. However, as n increases, the amplitude of oscillation decreases 
and the oscillation will eventually cease as T gets smaller and smaller 
and n gets larger and larger. 

If y and T grq are increased then the curves of -1/N(z) move up with 
the end point shifting according to Eq. (4-15). Figures 4-7 through 4-9 
show the plots of -1/N(z) for n = 2, 3 and 4, respectively, with y s 8470 
and T qf0 « 0.0424. In this case the lowest point is at -23.66 db and 


IR quality 
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Figure 4-7. Discrete describing function for the frictional nonlinearity, 
n » 2. 




















self- sustained oscillations are possible only for n ■ 2 over a fixed 
small range of sampling periods. If T Is chosen to exclude the critical 
region, self-sustained oscillations can be avoided. 
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5. Discrete Describing Function of A Quantizer 
5" 1 • Int roduction 

The effects of quantization and the study of quantization error 
conducted In Chapter 1 are all based on the assumption that the dlqltal 
system with quantization Is stable and free from sustained oscillations. 
Since a quantizer Is a nonlinear element, It can cause self-sustained 
oscillations. In Chapter 2 the condition of self-sustained oscillations 
In the analog low-cost LST system with quantizer Is studied by use of 
the continuous-data describing function. However, in reality, the low- 
cost LST system Is a digital system. The Interaction between the sampling 
operation and the quantizer will bring about phenomena which can be 
grossly different from that In an analog system. Therefore, It is 
essential that the discrete describing function (DDF) of a quantizer be 
derived. To the authors' knowledge the DDF of a quantizer has not been 
derived In the past. 

Figure 5-1 shows the Input-output characteristics of a quantizer. 

The Input of the quantizer is denoted by e*(t), and the output by y*(t). 

It is assumed that the input of the quantizer is the output of an ideal 
sampler. Therefore, e*(t) and y*(t) are trains of Impulses. Furthermore, 
it is assumed that the input to the sampler, e(t), Is a cosine wave, and 
thus the amplitude of e*(t) Is modulated by a cosine wave; that Is, 

e(t) * E cos (wt ♦ <J>) (5-1 ) 

00 

e*(t) * E l cos (kwT ♦ <J.)<S(t - kT) (5-2) 

k=o 





where E Is a constant, m Is the frequency In radians per second, T Is 
the sampling period in seconds, and 6 (t) denotes the unit Impulse 
function. The z-transform of e*(t) is 
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Eli) ■ -5 [ (z - coswT) cos<(> - sin wT sir $] 

Z - 2 Z COS cdT ♦ 1 

The z-transform of y*(t) is denoted by Y(z). 

The discrete describing function of the quantizer nonlinearity is 
defined as 



As in the case of the relay-type nonlinearity, we assume that 
because of the periodic nature of the sampler, e*(t) and y*(t) are all 
periodic functions of period nT, where n is a positive integer greater 
than or equal to two. Thus, 


or 


<aT 


2tt 

n 


n = 2, 3, 4, ... 


The DDF N ( z ) is incorporated in the "characteristic equation" of 
the system, 

1 + N(z)G eq (z) = 0 


(5-3) 


(5-4) 


(5-5) 


(5-6) 


(5-7) 


4 ‘ 

l 
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for the determination of the condition of self-sustained oscillations, 
where G 'z) denotes the linear transfer function which the quantizer 
nonlinearity sees. Graphical'*./, the condition of self-sustained 
oscillation characterized by the period nT is determined by the Inter- 
sections of the G e( ^(z) trajectories with the critical regions of -1/N(z), 
all for the same n. Therefore, the DDF problem Involves the determination 
of the critical regions of -1/N(z) for the quantizer for n » 2, 3 




5-2. The DDF of A Quantizer for n 2 


In order to Illustrate the derivation of the DDF of the quantizer 
nonlinearity, we shall first consider the case of n ■ 2; that Is, the 
self-sustained oscillation Is characterized by the period which Is equal 
to twice the sampling period. 

For this mode of oscillation, the waveforms of e(t), e*(t), and 
y*(t) are shown In Fig. 5-2. 

For n * 2, Eq. (5-6) gives u>T * n; Eq. (5-3) becomes 

E(z) • (5-8) 

The z-transform of y*(t) is written 

Y(z) “7TT (5 ' 9) 

where k is a positive Integer. In this case It is assumed that the 
value of E is constrained by the following equation: 


Using Eqs. (5-8) and (5-9), we have 

1 _ E cos <{> 

“ NTH kh 


(5-10) 


(5-11) 


For a given set of values of k, h, and 4>. the constraints on the 




values of E for n = 2 are 
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E - ( 2 , k - -t.Uil 
max 2 cos <t> 


E , . i ff . -■ Jin 

min 2 cos <t> 


Substituting E max and E m ^ n Into Eq. (5-11), we have 


m 




N{TT 


(k) - - V 


The last two equations Indicate that the critical region fo*' -1/N(z) 
for n * 2 is the line which extends from -(2k - l)/2k to -(2k + I )/2k 
on the negative real axis in the polar coordinates. 

It can be shown that 


N(W (, ' f1)< N(zT (k) 

7 maw ' 7 in-4 v 


Since 


N(^T 


(k+1 ) 


2 ( k + 1) + 1 

2(ic+T) 


2(k + IT 


2k + 1 
2k 


1 +2* 


thus, Eq. (5-16) is verified. 




Figure 5-3. Critical region of -1/N (z) for qua 
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Similarly, using Eq. (5-15) we can show that, in general, 




min 


(k+1) 


> NTiT 


ini n 


(k) 


(5-19) 


Thus, let S(k) be the set which is bounded by 


NTzT 


(k) 


and 


max 


NTzT 


(k) 


min 


and S(k+1) be the set which is bounded by 


N(zT 


(k+1) and 


max 


N(D 


(k+1) 


min 


then 


S ( k+1 ) C S(k) 


(5-20) 


for k= 1,2, 3, ... 

Equations (5-14) and (5-15) also show that 


tim 

k-*o 


m 


(k) 


max 


ton 

k-*o 


m 


min 


(k) 


= -1 


(5-21) 


This result implies that as the number of quantization levels increases, 
the quantizer characteristics approach a linear gain, and the critical 
region reduces to the (-1, jO) point in the complex plane. 

The critical region for n = 2 is shown in Fig. 5-3. 
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5-3. The DDF of A Quantizer for n = 4 

Before embarking on the derivations of general expressions of 
-1/N(z) for all n > 2, the case for n * 4 will be considered as an 
Illustrative example. 

For r. * 4, the input pulse train of the quantizer can assume a 
maximum of two different pulse amplitudes, k Q and k^. The mode of 
oscillation is thus characterized by A = (k Q , k^), where k Q and kj are 
positive integers. Figure 5-4 illustrates the input and the output of 
the sampler when the former is shifted through 360 degrees. It is 
observed from these waveforms, as well as will be verified later by 
equations, that the critical region for n * 4 repeats every 90 degrees 
for <J>, so that only the range for -45° < 4> < 45° needs be considered. 

For n * 4, cos wT ■ cos 90° = 0, z = j. Thus, Eq. (5-3) becomes 

E(z) " — (z cos ^ - sin4>) (5-22) 

z'* + 1 

With reference to Fig. 5-4, the expressions of the input to the 
quantizer for th’ various ranges of <f> are written as follows: 


h(k z + k,)z 

-90° < * < 0° Y(z) = ^ — (5-23) 

z L + 1 

h(-k z + k,)z 

-180° < d> < -90° Y(z) * % — (5-24) 

: ♦ 1 

h(k z - k.)z 

0° < $ < 90° Y(z) = 1 (5-25) 

z' + 1 


*» 


V 





0* < ♦ < 45° (k, < k 0 ) 



♦ - 45° (k, • k o ) 



45° < * < 90° (k, > k Q ) 


* • 90° (k o ■ 0) 


90° < ♦ < 135° (i. 0 < k,) 


♦ • 135° (k Q • kj) 


135° < $ < 180° (k Q > k,) 


« « 180° (k, • 0) 


Figure 5-4. Input and output signal waveforms of sampler preceding 
the quantizer as varies through 360 8 . 
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♦ • 0* (k, • 0) 


t -45* < ♦ < 0' (k Q > k,) 


♦ • *45* U 0 ' k,> 


-90' < 4 < -45' (k Q < k j ) 


♦ • -90° (k o • 0) 


-135' < * < -90' (k. > k ) 
■ o 


♦ • -135° (k, • k Q ) 


-180° < « < -135' (k, < k Q ) 


Figure 5-4. Input and output signal waveform of sampler preceding 
the quantizer as $ varies through 360®. 
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90° < ♦ < 180 


h(-k z - k,)z 

*<z) • % — 

z £ ♦ 1 


(5-26) 


The constraints on the magnitude of E over these ranges of i are 
tabulated as follows: 


(2k - l)h (2k ♦ l)h 

9—x < I E COS d> | < 


(5-27) 


and 


(2k. - 1 )h (2k, + 1 )h 

g < |E sin <f> | < ^ 


(5-28) 


These inequality conditions lead to maximum and minimum bounds on E. 

When tnese E„, w and i. are substituted into 
max min 


m m 




(5-29) 


and along with Eqs. (5-23) through (5-26), the equations that define 
the boundaries of the critical regions are obtained as 


1 


m 


and 


max 


NTzl 


min 


For instance, for -90° < <}> < 0°, from Eq. (5-27), 


(2k o - l)h 
'min 1 2 cos <J> ~ 


(5-30) 


Using Eqs. (5-22), (5-23), and (5-29), we have 
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nTzT 


2k o - 1 


■ oY'ir’ry TT"\ ( tan $ ■ J) 

mini V 


(5-31) 


Rationalizing the last equation, we have 


N(W 


2k - 1 

• 9 jr - [ (k, tan <J> - k ) + j(-k tan <* - k.)] 

mini 2(k 2 ♦ k 2 ) 1 0 0 1 


(5-32) 


Let 




min 1 


- R ♦ jl 


(5-33) 


where 


2k - 1 

— 9 — ( k i tan * ‘ k n ^ 

2(k 2 ♦ k 2 ) 1 0 


(5-34) 


2k ' 1 

1 “ 5 T (' k n tan< * ' k l) 

2<k‘ * k‘) 0 1 


(5-35) 


Then, solving for tan 4> from Eg. (5-34) gives 


R + k C 

Mn * • T[r~ 


(5-36) 


where 


?t o- 1 

2( k o * k f> 


(5-37) 
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Substitution of Eq. (5-36) Into Eq. (5-35) yields after simplification. 


k„ 2k - 1 

1 ’ • *7 R • -fcp 


which Is the equation of a straight line In the complex -1/N(z) domain. 
Similarly, from Eq. (5-27), 


max 1 


(2k o ♦ l)h 


cos <t> 


-90° < <j> < O' 


Following through the same procedure as described above, we have 


■2i%rrb ,un *- J> 


which is represented by a straight line with the equation 


2k_ + 1 


I - - r 2 - R - 


The boundaries for -1/N(z) when the constraint equation of Eq. (5-28) 
Is used are denoted as 


NtiT 


max 2 


1 

NUT 


min 2 


The equations which define the boundaries of -1/N(z) for the entire 
range of <j> are tabulated below: 


-90° < <t> < 0° 
90° < 4> < 1 80 c 


K 2k . - i 

_0 n o 

k ] K ‘ 2k 1 


(min 1 ) 


i 
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k n 2 K + 1 

I ■ • R * — (max 1 ) 



(min 2) 


2k. ♦ 1 

R + -k— 
0 


(max 2) 


(5-42) 


180° < (fr < -90° 
0° £ <J> £ 90° 



(min 1 ) 


(max 1 ) 



(min 2) 


k. 2k. + 1 

1 * ' r R ' — (max 2) (5-43) 

o o 

These equations show that for the oscillatory mode of n = 4 and 

A * (k Q , k^), the critical regions are bounded by straight lines, and 

that it is sufficient to consider the range of -90° < 4> <_ 90°. 

However, Fig. 5-4 shows that k > k, for -45° < a < 45° and k < k. 

o—i — — o—l 

for 45° < |<J>| £ 90°, thus, the critical regions for modes A * (k Q , k^) 
and A s (k^, k Q ) are identical. It is necessary only to consider the 
range of <J> from -45° to +45°. In general, one does not have to use the 
constraints on <j>, as the intersects of the eight equations in Eqs. (5-42) 
and (5-43) will naturally define the critical region. 
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Mode* A • |fe o , 0) and A ■ (0, fcj) 

Figure 5-4 shows that when ♦ ■ 0° or 180®, k^ ■ 0, and the mode Is 
described as A « (k Q , 0). Similarly, when ♦ ■ ♦90°, the mode Is 
A ■ (0. kj ). 

For the mode A ■ (k Q , 0), Eqs. (5-42) and (5-43) are reduced to 


2k 0- ’ 


(min 1 ) 


2k ♦ 1 
o 


(max 1 ) 


1 ■-«- 
o 


(min 2) 


I ■ jk- (max 2) I 

o 

These four lines define a square as shov/n in Fig. 5-5. 

For the mode A ■ (0, k^, It Is simple to show that the critical 
region Is described by equations of the same form of Eq. (5-44) with 
k Q replaced by kj . Therefore, the square region shown in Fig. 5-5 
is also for the mode A ■ (0, kj ) with k Q replaced by k^. 

Mode. A • (fe y , fej), k o • fej 

When $ * i45°, *135°, k Q ■ kj, as shown In Fig. 5-4. The critical 
region is now described by the following four equations: 


2k - 1 


I = -R 





I « -R - 


2k 0 ♦ 1 

"^T 


I 


* R ♦ 


2 k o- 1 

“*T 


I « R + 


2k + 1 

o 

2k o 


The critical region is again a square, as shown in Fig. 5-5. 

In general, the critical region for the mode A = (k Q , k^ ) is 
defined by the intersects of the eight equations in Eqs. (5-42) and 
(5-43). As an illustration. Fig. 5-6 show., toe critical region for 
A * (2, 1) which is also for A * (1, 2). Notice that this critical 
region is bounded by that of A * (2, 0) or A = (0, 2). 


(5-45) 
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(1) I 

( 2 ) 1 

( 3 ) I 

( 4 ) I 


-2R - 3/2 
-2R - 5/2 
2R + 3/2 
2R ♦ 5/2 


(5) I ■ R/2 ♦ 1/4 

(6) I - R/2 + 3/4 

(7) I » -R/2 - 1/4 

(8) 1 « -R/2 - 3/4 


\ / 

\ / 

\ N ' i 

\ A » (2,0) ^ / A * (2,2)/ 

-\A - (0,2) V / J' 

\ i ' N / 

\ 1 / 

• . . ' V - - 


-1 A\, 


✓ 


A * ( 2 , 1 /) 
A = ( 1,2 y' 


V 

/\ 

/ \ \ 


Figure 5-6. Critical regions of A = (2,0), (0,2), (2,1), 
(1 ,2) , and (2,2) , for n = 4. 
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